THE EULER-MARUYAMA APPROXIMATIONS FOR THE CEV 

MODEL. 
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_ , Abstract. The CEV model is given by the stochastic differential equation 

^^ . Xt = Xo + /q ^iXsds + /g a{X+)PdWa, 5 < P < 1. It features a non-Lipschitz 

diffusion coefficient and gets absorbed at zero with a positive probability. We 

^^' show the weak convergence of Euler-Maruyama approximations X" to the pro- 

j^ , cess Xt, < i < T, in the Skorokhod metric. We give a new approximation 

by continuous processes which allows to relax some technical conditions in the 
proof of weak convergence in I19| done in terms of discrete time martingale 
problem. We calculate ruin probabilities as an example of such approxima- 
tion. We establish that the ruin probability evaluated by simulations is not 
guaranteed to converge to the theoretical one, because the point zero is a dis- 
^^ ' continuity point of the limiting distribution. To establish such convergence we 

use the Levy metric, and also confirm the convergence numerically. Although 
the result is given for the specific model, our method works in a more general 

r^ ■ case of non-Lipschitz diffusion with absorbtion. 
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1. Introduction and the Main Result. 

We consider the Constant Elasticity of Variance (CEV) model (e.g. [3], [5]) 
defined by the fto equation with respect to a Brownian motion Wt 

Xt = Xo+ [ fiX,ds+ [ a{X+)PdWs, (f.f) 

Jo Jo 

f^ ' where a;"*" = OVx, constants p G [i f), yu £ M, cr > 0, and a positive initial condition 

Xq. Denote by r the first time of hitting zero r = inf{t : Xt = 0}. It is known that 
f~^ ' zero is an absorbing state, e.g. [TS], and that P(t < 00) > 0. 

f^ I The aim of this paper is to prove validity of approximations of expectations of 

some functionals of these diffusions by the Euler-Maruyama scheme. For example, 
in finance this model represents a price process and it is important to evaluate ex- 
pectations of payoffs depending on past prices E(y(Xro,T]); as well as for evaluation of 
. ^ the ruin probability P(r < T) by simulations. In population modeling the diffusion 

jrt ' with p = 5 represents the size of a population and is known as Feller's branching 

diffusion. In this case r represents the time to extinction of the population. 

When exact theoretical expressions of such functionals exist then a comparison 
of simulated and exact results show how good the approximation is. However, in 
many cases exact expressions are not available. Then one needs to justify such 
approximations. Weak convergence established below assures convergence of ex- 
pectations of bounded and continuous functionals of simulated values to the exact 
ones. 

To allow a larger family of functionals to be approximated we show the weak con- 
vergence of approximations in the Skorokhod metric rather than uniform. We also 
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consider approximation of the ruin probability P{t < t), which is the expectation 
of the past dependent and discontinuous in this metric functional I[r<T}- 

The fact that the diffusion coefHcient is singular and non-Lipschitz makes the 
analysis non-standard. Feller in |^ showed for the case p = 1/2 that the solution of 
the Fokkcr-Plank equation exists and is unique and gave its fundamental solution. 
This fact is used for evaluation of ruin probability by claiming that for any T, Xt 
has a density on (0, oo), and in particular its distribution function is continuous at 
any positive point. 

Previous various results on the Euler-Maruyama algorithm for different type 
of diffusion models can be found in Kloeden and Flatten [T!] and Milstein and 
Tretyakov [16] . They are concerned mostly with the approximation of the final 
value Xt, rather than the whole trajectory Xt, < t < T. When both the drift 
and diffusion coefficients are Lipschitz and the diffusion coefficient is nonsingular, 
the standard theory applies. There is a body of literature on the topic of approx- 
imations, see e.g. Bally and Talay [T], Bossy and Diop [3], Gyongy [7], Gyongy 
and Krylov [5], Halidias and Kloeden [3], Higham, Mao and Stuart [TI], Hutzen- 
thaler and Jentzen [T^] . Zahle [12] considers the case when the classical "Lipschitz- 
Lipschitz" setting fails. However, as a rule, the cases considered in the literature 
are such that the drift and diffusion coefficients exclude the absorbtion effect even 
when the diffusion coefficient is singular such as in Bcssel diffusion. 

We chose to give the results for the particular model of CEV for the sake of 
transparency, and because this model is of significance in applications. The reader 
will note that the proofs are sufficiently involved already in this special case. How- 
ever, our method and result hold in a more general case of non-Lipschitz diffusion 
with absorbtion. 

Recall next the Euler-Maruyama scheme for the diffusion process Xt- Taking 
for simplicity equidistant partitions of [0,T], = ^q < t" < . . . < t"^_^ < t^ = 
T, t^ — t'^_i = — , it is defined by the following recursion 

^tl = ^0, 

XJi = X- + ^,-{X-. )+ + aiX-.+ fJ^i,,, (1.2) 

k k — 1 fi k — 1 k — 1 y 11 

xr = x^._^, ie[i^l,i^), k = i,...n, 

where (^fc)fc>i is i.i.d. sequence of (0, 1)-Gaussian random variables. 

The process XJi can be considered as a discrete time semimartingale (see, e.g., 
Zahle [19]) and X" as a continuous time semimartingale with discontinuous paths. 
It may seem that there is not much difference between discrete or continuous time 
approximations. However, using a continuous time continuous approximation of X" 
introduced in the sequel enables us to simplify the proof of the weak convergence 
and avoid some technical assumptions used in [TO], e.g. (2.3) ibid. 

The paths of the process X" = (X")4g[o_T] are right continuous piece- wise con- 
stant functions with left limits belonging to the Skorokhod space B = D[o,t]- Con- 
sider a metric space (D,do) endowed with the Skorokhod metric do: if x, y G D, 
then 

do{x,y) ^ ini I sup \xt-y^{t)\+ sup '"" 



¥=6* I. tg[0,T] 0<s<t<l 



los 



t 



where $ is a set of strictly increasing continuous functions ip ~ i^it))o<t<T with 
ifiO) = 0, ipiT) = 1. 

Denote by Q" the probability measure on the space (D, i^), where & is the Borel 
(T-algebra, of the distribution of X" ~ {X")te[o,T], and by Q the distribution of 



X = {Xt)t^io^T]- Since X is a continuous process, Q(C) = 1, where C denotes the 
subspace of continuous functions. 

Evaluations of functionals by simulations is justified by the convergence of mea- 
sures in the Skorokhod space Q" — '^—^ Q, that is, 

n— )-oo 

lim / /(x)dQ" = / f{x)dQ 

for any bounded and continuous in the metric do function /(x). The aim of the 

paper is precisely to show this convergence. Note that such a function /(x) need 

not be continuous in the uniform metric g{x,y) = sup \xt —yt\- 

te[o,T] 

We remark on the use of the function x^ next. Since the process X = {Xt)te[o,T] 

is nonnegative, the notation Xj^ in (11. ip is used formally. However, it must be used 

for the approximating process Xl^^ in p.2p . because Xf can (and does) become 

negative. Since X" has piece-wise constant paths the stopping time 

t" = mi{t < T : Xj" < 0} 

much more likely corresponds to the negative value rather than the zero value of 
X" . The plus notation X""*" in p.2|) enables us to exclude negative values of X" 
everywhere except of X" . Thus, in simulations of f{X) we must use /(X"+). 
Next, the function g(x) := I{x^) inherits the properties of /(x) and is bounded 
and continuous in the metric do, moreover, when applied to simulations it converges 
to the desired limit 

lim E/(X"+) = lim Eg(X") = £g{X) = E/(X+) = E/(X). 

We introduce now the continuous approximation X" used in the proof of weak 
convergence 

X-^Xo + Y. fiiX-. )+ds + ^ / <j(X-+ rdWs, 



where W* is a Brownian motion such that Wt" — Wt" = \/—S.k- 

Denote Q" the distribution of (X")tgro ti- The main result is that Q" — - — > Q. 

Theorem 1.1. For any T > 0, the Euler-Maruyama approximation for the model 
p.2p converges weakly in the Skorokhod metric do to the limit process (Xt)tg[o.T] 
defined in p.l|) . 

The proof is done in three steps: first we show that the distance in the uniform 
metric between the Euler-Maruyama and the continuous scmimartingale approx- 
imations converges to zero in probability, secondly we show that the continuous 
approximation converges in the Skorokhod metric to the solution, and finally we 
deduce convergence of the Euler-Maruyama approximation. Thus the three steps 
in the proof are 

step 1. e(X",X")^^^0 



step 2. Q" -^^^ Q 

step 3. ^ , "^°° } => Q" -^^ Q. 

Q" " > Q ' "-'■^ 



(1.3) 



?l— >-oo 
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The proof of "step 1" uses standard stochastic calculus calculations for semimartin- 
gales. The proof of "step 2" follows from the results on diffusion approximation 
for semimartingales (see [15j Ch.8]). The proof of "step 3" is done as a standard 
application of Billingsley [2j Theorem 4.4, Ch.l, §4]. 

The paper is organized as follows. In Section 2 wc give preliminary results used 
in the proof, with some being of independent interest. The proof of the Theorem 
is given in Section 3. Section 4 gives simulations for the ruin probability P(r < T). 
To be self contained, the existence and uniqueness of solutions of the stochastic 
differential equation for the CEV model (|1.1|) is ahown in Section [5] 

2. Preliminaries 

We shall use the Doob maximal inequality for martingales (see e.g. [15[ Ch.l, 
§9] and ^] p. 201]) in the following context. Let (afe)o<A:<n-i and (^fc)i<fe<„ be 
random variables such that 

- Eal < oo, < fc < n- 1; 

- (6)i<fc<« is i.i.d. with E^i = 0, ES,l = 1; 

- ^fe and {ao, ■ ■ ■ , afe-i} are independent for any 1 < k < n. 

[nt] 

Set Mt = J2 (^j-i^j; where t G [0, 1] and [nt] = A: if < e (^, ^]. The process Mt 

0=} 

is a square integrable martingale for a suitable filtration. So, the Doob maximal 
inequality E(| sup^grg ii MA ) < 4EM^ is equivalent to 

k 2 '^ 

E.^i.^i^ |E"^-i^^ <4^Ea2_,. (2.1) 

Next we give the result useful in the analysis of products of random variables. 

Lemma 2.1. Let Xn, Yn be sequences of random variables such that X„ converges 
to zero in probability and expectations of Yn are uniformly bounded, sup„ E|y„| = 
r < cxD. Then XnYn converges to zero in probability. 

Proof. 

p(|x„r„| > e) - p(|r„| < c, |x„r„| > e) + p(|y„| > c, |x„y„| > e) 

<P{C\Xn\>e) + P{\Yn\>C). 
By the Chcbyshcv inequality P(|l^n| > C") < {j- Hence 



lim P(|X„y„|>e) <- ^0. 



n 



Next we need the following result of convergence to zero the expectation of 
double-maximum the Brownian motion increments over partitions. 

Lemma 2.2. Let Wt be a Brownian motion and {ii}, i ~ 1,. . . ,n an equidistant 
partition of [0,1]. Let W* = supjgrj._^ j.i \Wt ~ Wti_^\, and Af„ ~ maxi<i<„ W^j*. 
Then AIn converges to zero in h'' for any k > 1. 

Proof. We calculate fc-th moment of A/„, A: > 1. Since M„ > 

EAf,^ = / P(A/^' >x)dx^ / P(Af„ > x^^'')dx 
Jo Jo 

<n I P{wi>x^'^\dx. 



The inequality is obtained since {W*)i<i<n are i.i.d. random variables, therefore 
for any x > 



P(M„ > x^/'') = P( max W* > x^^'') = p( \J{W* > x^^''}) < nP(w^ > x^ 



Ik 



i=\ 



Next since sup^grg , /„, \Wt\ < sup^grg w„, Wt + sup^grg w„,(-Wt), it follows 



P(W* > x^/'') <P( sup Wt+ sup {-Wt) > x^/'^ 
<p( sup Wt> ia;!/^') +p( sup {-Wt) > ^x'^'') = 2pfiyi/„ > ^x^/^ 



sup ^ ., _ , 

.te[o,i/'i] ^ / Vte[o,i/?i]' ' 2 y V 

where we have used the well-known law of the maximum of Brownian motion. 
Thus 

EM„^ < 4nJ°°p(w,^^ > ^x'/Adx < in J°° p(\W,/,-,\ > ^x'/^dx 

= 4n2''E\W,fX = Cku^-'^/^ 

where Ck depends only on k, {Ck = 2'=+2E|^|'= with ^ ~ N{0, 1)). Hence il/„ 
converges to zero in L*^ for any fc > 2. But since convergence in L^ for a p > 1 
implies convergence in L*"' for any k E [l,p], the statement is proved. D 

The next elementary inequality is new and is instrumental in the proof. 

Lemma 2.3. For any x,y gR, and p G [51 l)j 



(x+)2P - {y+)^P\ < 



{2+\x\ + \y\)\x-yY>, pe(i,l). 



Proof. For p = i we have 

k"^ -y^\ = \x- y\I{x>0,y>0} + \x\hx>t),y<0} + W {x<t),y>0} 

<\x- y\I{x>o,y>Q} + \x- y\I{x>o.,y<Q} + \y - x\hx<t),y>a} < \x - y\. 
For p £ (i, 1) we have 



(2.2) 



|(^+)2P _ (y+)2p| ^ |(^+)p _ (y+n|(x+f + (y+)P| 
<(2 + N + M)|(x+r-(y+f|. 

Next we now show that 

\x^ -y^\ < \x~y\^- 

For X = y = it is obvious. 

Consider x > and y < 0. Then, taking into account the proved inequality 
|x+ — J/'*' I < |x — 2/1 (the statement of this lemma for p — ^) and the fact that 
y+ =0, we obtain 

|(x+r - (y+n = \x+\P = |x+ -y+\P<\x- y\P. 

Clearly, the inequality remains true for x < and y > 0. 
If both X and y are positive and x > y, then 



xP-yP\^/P = x 



^_.y_^p 



i/p 



< X 



l-(^ 



1/p 



< X 



1- 



(!) 



= \x~y\ 



and, in turn, |x^ — yP\ < |x — y\P. It is easy to see the inequality remains true for 
y > X. D 
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3. Proof of Theorem ll.il 

In the following result we show that the maximum of discrete approximations 
have uniformly bounded seconds moments. Henceforth r denotes a generic positive 
constant independent of n with different values at different appearances. 



Lemma 3.1. Let X^n he the Euler-Maruyama approximation defined in (|1.2[) . 



Then E max \XIi^\'^ < r. 

l<fc<Tl fc 



^n 12 

Proof. First wc bound from above max E|X^p. By (|1.2p 



E\X-.\^ = E\X?. |2(i + ^)\^2e(^„+ )2pZ:. (3.1) 

k k-l V 7j / fc-l 7j 

We use repeatedly the following bound of {X"^ )^p. Since (a;+)^P < |xp^ and 
2p < 2, we have (x+)^p < 1 + x^. Hence 

EiX^l/^<l + E\Xllf. (3.2) 

Next, for sufficiently large n, there exists r such that 

(1 + — )'<! + -. (3.3) 

n ' n 

Hence for sufficiently large n, we obtain from (|3.ip by using p.3p the recurrent 
inequality 

E|xr,p<E|xiJ._j^(i + i) + ^(i + E|xr,p 

= E\Xli \^(l+'-±^)+^:^E\X?. Kl + ^ 

*=-! V n J n fc-i V n/ 71 

Iterating it, for k < n wc obtain 

E|x5P..j(..i)Vit(i.£'-^" 



< 

Thus 



1 + -) =(X2 + r)0(e-^). 



max E|Xj"„p <r. (3.4) 

l<k<n I' 



From the definition of the scheme (|1.2p we obtain by iterations 



^ '^ / rp k 



max IX.'^I < |Xo| + ImI-E 1^"" J + \/- ,5^f.^ E'^(^"^,)'^^- 



J=l 



J=l 



Next we use the Cauchy-Schwarz inequality (X]i;=i "^i)^ ^ 'X]i=i'^i with / = 3. 
Proceeding from above we have 



E max IX" P 

l<A;<n k 



<3|Xor + 3LE(^|X,'J. l) +3-Emax |Ea(X,'^+)PO • (3.5) 






The Cauchy-Schwarz inequality ( X]?=i l-^t" I ) — "-X^^^il-^t" )^ coupled with 
the bound on the largest second moment shown above (j3.4p implies that the second 



term in p.Sp is bounded by a constant r independent of n. The last term is bounded 
by using the Doob's maximal inequality (j2.1[) and the bound in p.2[) 

k ^ n n 

J=l J=l j=l 

Using (|3.4p we can now see that the last term in p.Sp is bounded by r. D 

We proceed now to prove the steps in p. 31) . 

3.1. Step 1. The distance between the two approximations in the sup norm con- 
verges to zero in probability and in L^. 

Lemma 3.2. Eg2(X",X") > 0. 

Proof. From the definitions of X" and X" it follows that at the points of the 
partitions both approximations coincide and at intermediate points the following 
holds 



^i{X-.J+ds+ aiX^lJdWs. te[i^i,t^]. (3.6) 

By the formula (|3.6p . 

sup \x- - X- I < i^ |xr„ ^ I + a sup I w^* - Wt^^ ^ I ixr„ ^ r. 

Consequently, 

g(X", 1") = sup |Xr - Xn = max sup \X^ - X" | 

I IT^ 

<■'—'— max |X"„ I + (T max sup \Wt — Wt" \ max |X, 



fc-i ' 



n l<k<n fc-i l<fc<ri(gr(n (,.i ' ''^^ l<k<n k-l' 

The first term converges to zero in L^ using Lemma [XT] as Emaxi<fe<„ \X"n \ < r 
For the second term use Holder's inequality with parameters - and 2Z~- 

E( max sup \Wt - Wf^ I max I X" P ' 

l<fc<ntg[t„_^_i„]l '=-i'l<fc<n' '-"-i' ; 

2 \ 2-p 



< E 



max sup \Wt — Wt^_ \ 



n |2 



E max IXJ 

l<fc<n fc-i 



By Lemma [3.11 sup^. Emaxi<fe<„ jX^" p < r. Since i^ > 1, by Lemma 



2/(2-p) 



lim E( max sup \Wt — Wt" |) 



3.2. Step 2. Weak convergence of continuous approximations. 
Lemma 3.3. Q" -^^ Q. 



n 
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Proof. The proof rests on a general result on the weak convergence of seniimartin- 
gales to a diffusion [T2|, Theorem 1, Ch. 8, §3. This theorem states that for weak 
convergence to a diffusion it is enough to check convergence of the drifts and qua- 
dratic variations evaluated at the pre-limit processes. The processes X and X" are 
semimartingales with following decompositions 



Xt — Xq 



[ pLXsds+ [ aiX+)PdWs 
Jo Jo 



--Bt{X 



■■=Mt{X) 

n tAt" 



X- =Xo + Y. f^X^l^ds + ^ / aiX^l^ rdW, 



=B^iX") 



:=A/,"(X") 



where wc have denoted above 



(B) Bt{X) and Sf (X") are drifts 

(M) Mt{X) and il/"(X") arc continuous martingales with predictable quadratic 
variations 

{MUX)= f a^Xtf^ds, 



(M"),(X") = ^/ \\X-lf^ds. 



The above mentioned Theorem 1 of [15] , adapted to the present setting, states that 
the weak convergence takes place if the following three conditions hold. 



(a) Equation (II. ip has a unique (at least weak) solution 



prob. 



>0 



(b) sup |Bt(X")-S,"(^") 
te[o,T] 

(c) sup |(M)(X")-(M")t(X")|^^i^^O, 
te[o,T] 



> ^Q 



prob. 

n— )-oo 



n do_^Q, 



Wc proceed to verify these conditions. The existence and uniqueness of \l.\-\ is 
known, e.g. [IHI, and is also given for completeness Proposition 15.11 Hence (a) 
holds. To show (b) write, taking into account (|2.2p . 



sup 

te[o,T] 



B,{X^) - B^{X-)] < ImI E / 1^" - ^*L, 1^* ^ T\t,\Q{X^, X-) 



fe=i'"^fc-i 



Hence (b) holds by applying Lemma IX^ 
To prove (c) write the bound 



sup 

te[o,T] 



{MUX-) - (Af")*(X") <a'Y. K^"^)'" - (^*r.)''l^* 

= <t'J2 liX^f ~ {X^lflds, (3.7) 



k=l -"^fe-i 



where we have used that the two approximations coincide on the grid. By applying 
Lemma 12.31 we have further bound on the expression under the integral 



ixr)2^-(x"+ ) 



2p 



fc-l 



< 



2+ sup |Xr|+ sup \Xl^\) sup |X,"-X" |P, pe(i,l) 
te[o,T] te[o,T] ^ se[t^_i,tj] 



Hence for p = i the bound in p.7p becomes 



and the statement foUows by Lemma \WI^ 
For p E (^,1) the bound in p.7p becomes 



< ^2(2 + sup |Xr| + sup \X^\) E / sup \X: - XlTds 

^ tG[0,T] te[0,T] ^ fc=i"'tL!_iS6[t2-i.tfcl 

<Ta2(2+ sup |X,"|+ sup |X,"|)eP(X",X"). 



te[o,T] te[o,T] 

By Lemma |3. II Lemma l3.2l we have the product of two terms, one of which has 
uniformly bounded expectations and the second converges in probabihty to zero. 
By Lemma |2. II the product converges in probabihty to zero. 

Thus the conditions of the Theorem 1 of [15j are verified and weak convergence 
is proved. D 

3.3. Step 3. Weak convergence of the Euler-Maruyama approximations. 

Lemma 3.4. For any bounded and continuous in the metric da function f{x), 

U^ |E/(X") - E/(X)| < li^ |E/(X") - E/(X)| = 0. 

n— >-cxD n—yoo 

Proof. The result follows from the triangular inequality 

|E/(X") - E/(X)| < E|/(X") - /(X«)| + |E/(X") - E/(X)|. (3.8) 

Since convergence in the uniform metric imphes convergence in the Skorokhod met- 
ric, by Lemma 13.21 we have, since / is continuous in this metric, 

hm E|/(X")-/(X")HO. 

Taking now the hm sup in p.8p and using the previous result of weak convergence 
of X" to X proves the step 3. D 

4. Evaluation of ruin probability on a finite time interval by 

simulations. 

In this section we evaluate numerically a ruin probability P(r < T) on a time 
interval [0,r], where r = {i : Xf = 0} by Euler-Maruyama approximations. An 
explicit formula to the ruin probability, as a function of arguments p, Xq, /^, cr, T, is 
known in the case p = -^jT = oo: 

P(r < oo) = exp r ^Xo 

(see, e.g. [T31 P- 354]), but not for finite T. 



10 
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Naturally, P(r < T) is important in applications, and we study it for different 
values of parameters p, Xq, fi,a,T: 



Xo 



j_ 

10 

1 

4 
1, 



1 



M = 



-1, 



(7= 1,T = 



We combine the Euler-Maruyama simulation algorithm and the Monte-Carlo 
technique with 10^ runs per point. 

The basis for analysis is an obvious formula P(r < T) = P{Xt = 0). It allows us 
to deal with the distribution function F{x) := P{Xt < x) of Xt instead of harder 
to compute distribution function of r, P(r < T). Notice that 



Fix) 



0, a; < 

F(0) == P{Xt = 0) > 0, a- = 0, 



that is, F{0) — F{0—) > and so the distribution function F{x) has an atom at 
the point {0}. 

The measure Q is supported on the space of continuous functions. So the weak 



convergence of processes and measures Q*^ 



do 



-> Q implies weak convergence of 



law 



finite dimensional distributions, and in particular marginals, X^ > Xt- That 

n— >cjo 

is if F„{x) ~ P(X^ < x) then lim„_j.oo Fn{x) — F{x) at any point of continuity of 
F. Unfortunately 0, our point of interest, is an atom of F and we can not claim 
that 

lim F„(0) =F(0) = P(r <T). 

n— >oo 

In view of this uncertainty, we give approximations for lower and upper bounds 
of F{Q) by using the Levy metric (see e.g. [TU]'): 

C{Fn,F) = mi{h > : F„(a; - h) - h < F{x) < F„(a; + h) + h; V x}. 

It is known that weak convergence of distributions implies convergence in the Levy 

metric Q" -^ Q ==> lim C{Fn, F) = 0. 

n— >oo 

X = and a small suitable e„ , in each case determined experimentally, such that 
the interval [Fn{—en) — £n, Fn{en)+Sn] is small enough and declare that an estimate 
of F{0) belongs to this interval. Such values of £„ are pointed out below as last 
values in each Table. 

Simulations for various values of parameters show good fit of this procedure. 



Example 4.1. p 



4, a = l, Xo 



T - 



e 


P(X? <-e)-e 


P(X? < e) + e 


P{Xt = 0) 


2-10"^ 


.9679 
.9700 

.9738 


.9738 
.9738 

.9738 


[.9679, .9738] 
[.9700, .9738] 

[.9738, .9738] 


10-« = e„ 



Table 1. 
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Example 4.2. p 



i t^ 



+ l,a = l,Xo^^„T^9 



e 


P(X? <-e)-e 


P(X^ <e) + e 


P{Xt = 0) 


2-10"" 
10"" 


.8166 
.8174 

.8182 


.8192 
.8192 

.8192 


[.8166, .8192] 
[.8174, .8192] 

[.8182, .8192] 


5-10-^ = e„ 



Table 2. 



Example 4.3. p = \, ii = +1, a = I, Xo ^ \, T = 



£ 


P(.Y? <-e)-e 


P(X? <E) + e 


P(Xt = 0) 


2-10"" 

io-« 


.5960 
.5964 

.5970 


.5970 
.5970 

.5970 


[.5960, .5970] 
[.5964, .5970] 

[.5970, .5970] 


5-10-^ = £„ 



Table 3. 



Example 4.4. p = |, ^ = +1, a = 1, Xq = 1, T = 9 



e 


P(X? <-£)-£ 


P(X? < £) + e 


P{Xt = 0) 


3-10-'' 
2-10"^ 


.1344 
.1346 

.1346 


.1348 
.1348 

.1348 


[.1344, .1348] 
[.1346, .1348] 

[.1346, .1348] 


lO-'' = e„ 



Table 4. 



Example 4.5. p 



M = +l, cr = l, Xo = i,T = 9 



c; 


P(X? <-£)-£ 


P(X? < e) + e 


P(Xt = 0) 


3-10-« 
5-10-^ 


.6040 
.6040 

.6180 


.6206 
.6126 

.6206 


[.6040, .6206] 
[.6126, .6206] 

[.6206, .6206] 


2.5 • 10-^ = e„ 



Table 
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Example 4.6. p = |, ^ = +1, a = 1, Xq = 3, T = 9 



e 


P(X? <-e)-e 


P(X? < e) + e 


P(Xt = 0) 


3-10"* 
5-10-9 


.3794 
.3816 

.3838 


.3864 
.3864 

.3864 


[.3794, .3864] 
[.3816, .3864] 

[.3838, .3864] 


2.5 - 10-9 = e-a 



Table 6. 



Example 4.7. p = |, /x = +1, a = 1, Xo = 1, T = 9 



e 


P(X? <-£)-£ 


P(X? < e) + e 


?{Xt = 0) 


2 -10-* 

io-» 


.0752 
.0768 

.0782 


.0790 
.0790 

.0790 


[.0752, .0790] 
[.0768, .0790] 

[.0782, .0790] 


2.5 - 10-9 = e„ 



Table 7. 



Example 4.8. p = |, ^ = -1, a = 1, Xo = §, T = 3 



£ 


P(X? <-£)-£ 


P(X^ < £) + £ 


P(Xt = 0) 


10-** 
5-10-9 


.8536 
.8700 

.8757 


.8803 
.8803 

.8803 


[.8536, .0790] 
[.8700, .0790] 

[.8757, .0790] 


2.5 - 10-9 = £„ 



Table 8. 



5. Existence and Uniqueness of solution in the CEV model 

Proposition 5.1. The Ito's equation (|l.ip possesses a unique strong nonnegative 
solution. 

Proof. Define a sequence of processes indexed by integers i, i > l/^o, {Xl)t>o, 
sucfi tliat XI is a strong solutions to tlie following sde 



dXl = nXldt + a{i-^ V \Xl\)PdWt, X'„ = Xa 



(5.1) 



The diffusion coefficient a{i ^ V |a;|)P is Lipschitz continuous therefore X^ is the 
unique strong solution of (|5.ip . Set -di = inf{i : X]: = i^^}. Note that for t <di 

Xl+^ = XI 

and it follows that i?i+i > ^,;. A strong solution to (jl.ip is constructed by a natural 
prolongation 

i>n 

Finally, Yamada-Watanabe's theorem (see, e.g. [T71 Rogers and Williams, p. 
265]) guaranties the uniqueness of the strong solution of equation (jl.ip . because 

D 



the Holder parameter p > 5 
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